Stabilizing dynamical localization in driven tunneling 
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Dynamical localization of a tunneling system by means of an oscillating external field is examined for an 
arbitrary doublet of tunneling states. We show that the condition of an exact crossing of Floquet levels, required 
for localization in a symmetric system, can be substantially relaxed for tunneling systems with broken symmetry. 
This generalization equally applies to tunneling systems coupled to a dissipative environment. 

PACS numbers: 03.65.-w, 05.30.-d, 61.16.Ch, 62.65.+k 
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A tunneling particle can be localized indefinitely by an ex- 
ternal periodic force. This surprising prediction was made re- 
cently by Grossmann et al. With properly chosen fre- 
quency and amplitude the external force can pin the tunneling 
particle, and the amplitude of the tunneling oscillations can be 
made arbitrarily small. Numerous methods have since been 
applied to study this effect, focusing on special cases which 
lend themselves to analytical solutions or on their gen- 
eralizations to tunneling systems with dissipation [p|~p"oj]. Pre- 
vious works are based on the central assumption that an exact 
crossing of Floquet levels is necessary for dynamical localiza- 
tion to occur. For Floquet levels to cross, the tunneling system 
typically has to be symmetric, and the driving amplitude must 
be selected from a set of discrete "magic numbers", different 
for each driving frequency |[l^] . 

In this Letter, we show that the tunneling particle can re- 
main pinned even if both of the above conditions are relaxed 
simultaneously. This implies that the class of systems in 
which the dynamical localization can potentially be detected 
is much wider than previously thought. These include tun- 
neling systems which are not intrinsically symmetric, such 
as tunneling of an atom between a sample and an STM tip 
[ 16 ], or tunneling in disordered systems like structural glasses 
[ 1S,^0|]. We demonstrate that these findings also apply to dis- 
sipative tunneling systems, where the localization effect mani- 
fests itself through a suppression of tunneling oscillations and 
a drastic slowing down of the incoherent dynamics. 

To study driven tunneling systems, we can either use a 
double-well potential model [|lj) or, more simply, a truncated 
two-state Hamiltonian [[|] 
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This truncation is valid if all parameters are small compared 
to the lowest oscillation frequency ui osc associated with the 
double well, and for a barrier height larger than ?iui osc . a x 
and a z are Pauli spin matrices, hA is the tunnel splitting, hso 
an intrinsic energy bias, and he sin u>t its modulation by an 
external periodic force. 

Despite the simplicity of this Hamiltonian, its dynamics 
shows surprising effects in parameter regimes for which the 
rotating-wave approximation fails. A general analytical so- 
lution is not available, hence a detailed understanding of its 
dynamics must rely on numerics. 



The Hamiltonian ([!]) is non-conservative, but it exhibits dis- 
crete time-translational invariance with period T = 2it/u). 
Consequently, the one-period time translation operator Ut = 
U(T, 0) contains all essential information about the dynam- 
ics. Previous studies [|l]-§] concentrated on the case where two 
Floquet states, the eigenstates of Ut in the two-state descrip- 
tion, are degenerate. In that case, Ut is the identity operator 
(or its negative), and the time dependence of all observables 
is strictly periodic with period T. Within the interval T, the 
population p(t) of the initially occupied site exhibits small os- 
cillations, with a minimum p m j n such that |1 — Pminl <^ 1- 

In this work, we will consider this inequality the definition 
of localization and make no assumptions about the spectrum 
of Ut- In general, Pmin is to be defined as the largest lower 
bound (infimum) of pit) for any positive t, 
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Although t is unbounded, p mm can be evaluated without un- 
controlled approximations or extrapolation. Using the group 
properties of discrete time translations, the dynamics at arbi- 
trary time t can be constructed from the time translation oper- 
ator at times t < T, 



u t = u T m 



(3) 



Here we have defined Ut — U(t, 0), with t = nT + t, 
< r < T. The lower bound p m i n is thus given by 



Pn 



min inf p n (t) 

0<T<Tn>0 



where 



p n [r) = \(MUtU%\^)Y 



(4) 
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generally defines an infinite set of functions on the interval 

< r < T. 

To obtain the infimum with respect to n, it is useful to rep- 
resent the SU(2) matrix Ut by a generator g and an angle ipo, 
Ut = exp(^(fog) . The unit vector n in g — n ■ a denotes 
the rotation axis of the corresponding SO(3) rotation. Now p n 
reads 



P„(t) = \ (ipo\U T eiq)Un(pog)\ip ) 



(6) 
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and one can define a quantity closely related to p n (r) with a 
continuous index <p, 

P V ( T ) = \(^o\U T exp(%(pg)\i/j )\ 2 , (7) 
where < p < Air. The set of functions {p n (T~)} is a subset 

min p v {r) < inf p n (r). (8) 

0<ip<47T r n>0 

For generic parameters, ipo /tt is an irrational number. The set 
{exp(|n(/3 s) | n > 0} is a dense subset of {exp(|</?g) | < 
<p < 47r}, and strict equality holds in But in the degen- 
erate case of rational ipo / tt, the situation is different. p(t) is 
periodic, and there are only a finite number of functions p n (r) 
that are distinct. We shall consider the latter case irrelevant in 
the long-time limit, because any minute drift or other aberra- 
tion in the driving frequency or amplitude will reduce it to the 
generic case. We will concentrate on the generic case in the 
following. 

For the two-state system, p v (r) is a continuous function of 
(p. and r, taking the form p v (r) = a T + b T sin ip + c T cos tp. 
The coefficients a T , b T , and c T are algebraic functions of the 
matrix elements of Ut and U T . From the numerical solution 
of the Schrodinger equation over a single period one can thus 
easily obtain the entire set of functions p v (r), and also the 
lower and upper envelopes of the set of curves p n (r), 

pie(r) = inf P„(r) = min p v ( T ), (9) 

n>0 0<<p<2tt 

Puc(t) = supp n (r) = max p v {r). (10) 

n>0 0<v<2tt 

According to (Q), p m i n is then given by the minimum of the 
lower envelope, 

Pmin = min p\ c {r). (11) 

0<r<T 

For a symmetric system driven at frequencies uj 3> A, p m i n 
has been shown to be almost unity if the ratio i/u> belongs to a 
set of discrete 'magic numbers' for which Ut is diagonal and 
its eigenstates degenerate OU^]. In the high-frequency limit, 
these numbers are the zeros z n of the Bessel function Jq(z). 
But with even a slight deviation from these numbers, the de- 
generacy of the Floquet states will be lifted and Ut will have 
delocalized eigenstates, which manifest themselves as slow 
oscillations in pit) with large amplitude and a reduced tunnel- 
ing frequency AJq(£/oj). 

The central finding of this Letter is that the instability of 
the localization effect can be suppressed by imposing a small 
static bias £o on the tunneling system. This is demonstrated by 
the numerical data in Fig. la. For a symmetric system (dashed 
curve), a slight detuning of the amplitude from its magic value 
leads to to a slow coherent transition that depopulates the ini- 
tial state. Adding a small bias, however, prevents this depopu- 
lation even out to very long times (solid curve). Fig. lb shows 
the corresponding set of curves p n {r) for < n < 3 as well 



as the lower and upper envelopes p\ c (t) and p uc (t). This rig- 
orous result for p\ c shows unequivocally that a very robust 
localization persists indefinitely for the detuned system. The 
periodic function p(t) for the corresponding symmetric sys- 
tem with perfectly tuned amplitude e/ui « 2.393 is given for 
comparison (dashed curve). Both cases show roughly equal 

Pmin- 

How large does the bias need to be in order to effect dynam- 
ical localization, and is there a maximum allowable value for 
it? The precise answer depends quite sensitively on the driv- 
ing parameters. We can look for the answer using the seminu- 
merical procedure outlined above. The global minimum in Eq. 
dTT| ) can easily be obtained by linear search because pi e (r) is 
smooth. Fig. 2 presents quantitative results about the depen- 
dence of p m i n on £o for different £ in a tunneling system driven 
at uj = 5A. The width of the minimum at Eq = vanishes as 
e/u) approaches the localization point of the symmetric sys- 
tem, i.e., the bias required to enable localization goes to zero. 
On the other hand, for very large bias £o ~ uj, p m i n drops 
drastically, and the dynamics is then characterized by large 
oscillations [S]. 

A qualitative understanding of this generalized localization 
effect is obtained by evaluating the dynamics at finite time 
analytically in the limit of high frequency [^|^]. For uj 3> 
max(eo, A), we find |i^o| "C 1, and g — N~ 1 (AJq(£/uj)<t x + 
sqcTz). From this we conclude that localization is possible if 
either (i) <po = (Floquet level crossing) or (ii) | Jo(£/u>)\ <C 
£o/A| (axis n aligned with the z axis). The two cases are 
closely related. In either case, driving modulates the transi- 
tions amplitudes giving rise to the off-diagonal elements of 
Ut- In the first case, this is done in such a way that destruc- 
tive interference makes the off-diagonal elements vanish pre- 
cisely. In the second case, this cancellation is imperfect, but 
the residual off-diagonal elements are small enough to permit 
localized eigenstates of Ut- 

In a condensed matter environment, tunneling is coupled 
to fluctuations of the surrounding medium, and the delicate 
cancellations underlying this localization effect are disturbed 
by dissipation and dephasing. The standard approach to lin- 
ear dissipation resulting from such a coupling generalizes 
the two-state system to the driven spin-boson Hamiltonian 

IB, 

H = H a + ^u>„(ala v + |) + y ^C„(a„ +at)a z . 

Here and in the following we set h = 1; qo is the distance 
between tunneling sites. The effect of the harmonic environ- 
ment is fully characterized by the spectral density J(ui') — 
n Cl ^(y' ~ Wi/)- Its specific form depends on the den- 
sity of states of the dissipative environment and the details of 
their local interaction with the tunneling system. The Ohmic 
form J(oj') — rjuj' exp(— uj'/u) c ) is one of the most frequently 
investigated [ pTp2[ ]. qo and r/ appear in the dynamics only 
through the dimensionless dissipation constant a = riq^/Zir. 
In the case of large uj c applicable to tunneling in solid-state 
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systems, A and u c enter only through the scaled tunneling 
frequency A r = k(/\/uj c ) a/(1 ~ a) ■ 

With the notable exception of the noninteracting-blip ap- 
proximation (NIBA) [ pTp^ ] and certain limiting cases, the 
dynamics governed by the driven spin-boson Hamiltonian ap- 
pears to be intractable without resorting, at least partially, to 
numerical methods. The recently introduced chromostochas- 



tic quantum dynamics (CSQD) algorithm [23 1 makes possible 



the exact and efficient computation of all dynamical quanti- 
ties of interest in the spin-boson system. The CSQD method 
minimizes memory effects inherent to quantum dissipation by 
direct sampling of the colored noise produced by the environ- 
ment. The statistical error increases only moderately with t; 
in the examples given below it never exceeds 7 x 10~ 3 . 

The dynamics at low temperature and weak damping 
closely resembles that of the undamped system. Fig. 3 illus- 
trates dynamical localization in a weakly damped asymmetric 
tunneling system with e/uj — 2.3 driven at u> = 5A r . For 
a symmetric system, the same parameters result in a coher- 
ent population transfer between the two localized states. This 
transfer is strongly suppressed when a moderate bias So = — 1 
is introduced. The system now exhibits a gradual population 
decay on a much longer timescale. Comparison with the free 
population decay of the biased system (without driving) shows 
that the bias £o alone is not sufficient to cause localization. In 
this example, the tunneling particle is actually localized on the 
site whose energy is raised by the bias. Fig. 4 shows the same 
driven system with £o = — 1 (top) and eq = 1 (bottom) with 
little difference between the two curves, reminiscent of the 
symmetry seen in Fig. 2. The center plot shows a symmetric 
system driven at the same frequency, but with the amplitude 
needed for a Floquet crossing, if to w 2.393. Evidently, the 
two cases of localization discussed above lead to nearly iden- 
tical behavior also in the dissipative case. 

We have demonstrated that the remarkable effects resulting 
from dynamical localization in tunneling systems occur for a 
wide range of parameters that do not obey the restrictive con- 
ditions for an exact Floquet level crossing. We believe that 
these findings should stimulate further theoretical and exper- 
imental investigation. The prospect of an experimental real- 
ization has improved appreciably. The hard-to-satisfy require- 
ments of an extremely precise control of the driving ampli- 
tude and its perfect homogeneity throughout a sample can be 
dropped. The range of possible applications is substantially 
widened to include asymmetric and even disordered systems 
as candidates. Considering, e.g., the density and sound veloc- 
ity of vitreous silica and the strain coupling 7 2eV of its 
intrinsic tunneling systems [fZ4[], dynamical localization may 
be achieved by ultrasonic driving at low temperature. At a 
frequency of 1 GHz, the required amplitude for localization 
will be reached at an acoustic intensity of about 4 mW/cm 2 . 
Another phenomenon recently predicted and proposed as a 
testbed for theories of dissipative tunneling is the tunneling 
of a single atom between an STM or AFM tip and a sample 
fll8|]. The adsorbed atom forms dipoles of opposite direction 
on the sample and tip surfaces and thus experiences a force 



from an applied DC or AC tip voltage. For the Ni:Xe tunnel- 
ing system proposed in [ |l8| ] and characterized in more detail 
in [fZ5|], the localization condition for a 10 GHz driving fre- 
quency will be reached with an AC tip voltage of only about 
2 mV. 
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FIG. 1. (a) Time-dependent population p(t) for symmetric 
(dashed line) and asymmetric (eo = A, solid line) tunneling systems 
with u = 5A, and i/u> = 2.38. (b) Asymmetric system (solid lines): 
Non-periodic trajectories p„ (t) (thin) with periodic upper and lower 
envelope (thick) vs. At. Symmetric system (dashed line): Periodic 
p(t) of the dynamically localized system at i/u> ~ 2.393. 
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FIG. 2. Lower population bound p m i n vs. bias eo for various ra- 
tios ejus at u — 5 (unit A = 1). 
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FIG. 3. Population decay for weak dissipation (a — 0.01, 
zero temperature) with driving at ui — 5A r , e/co = 2.3, and 
u>c = 200 A r . Top curve - biased system: localized, long population 
decay time. Middle curve - symmetric system: coherent population 
transfer. Lower curve - biased system, no driving: fast tunneling 
oscillations and incoherent decay. 




FIG. 4. Slow population decay of a localized tunneling system 
driven at u — 5A r with weak dissipation (a = 0.01, zero tempera- 
ture). Top: eo = — 1, s/w = 2.3. Center: eo = 0, e/ui ~ 2.393. 
Bottom: eo = 1, i/u) ~ 2.3. 
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